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Abstract 

We present an analysis for the measurement of the Higgs Boson 
branching ratio BR(H — )• cc) for a hght Standard Model-Hke Higgs 
boson produced at 250 GeV centre of mass energy at the International 
Linear Collider (ILC). The tools and technique used for the analysis 
are described and relative statistical uncertainities are determined. 

1 Introduction 

The measurement of the Higgs absolute branching ratios to all possible 
species is an important part of the ILC program, giving a precision test 
of the Standard Model prediction that the Higgs boson couples to each par- 
ticle in proportion to its mass [Ij. These measurements also discriminate 
between different 'Beyond the SM' scenarios [21 13]. 

This study considers the Higgs decay mode to charm quarks which results 
in two- and four-jet final states. The expected SM Higgs boson branching 
ratio to charm quarks is equal to 3.09% for rriH = 120 GeV and one of the 
main difficulties of the analysis is to separate this signal from the background 
of Higgs decays to b-quarks which has substantially larger Br of 67.92%. 
The selected decay modes exercise the tagging of charm quarks which is 
particularly sensitive to the vertex detector performance 

2 Event Generation and Detector Simulation 

The signal events comprise of Higgs bosons decaying into charm quarks pro- 
duced in Higgs-strahlung, e+e~— )■ ZH ignoring other contributing processes. 
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The CP-even Higgs can also be produced via WW(ZZ)-fusion processes, 
e'^e~^uuW*W*—^uuH{e^e~^e^e^Z*Z*^e^e^E) as shown in Figure [H 
At low centre of mass energies the Higgs-strahlung process is dominant if 




Figure 1: Higgs production in (a) Higgs-strahlung, (b) WW-fusion and (c) 
ZZ-fusion processes 

Mh < 2Mw due to phase space considerations. The production cross-section 
is shown in Figure [2] [5]. At larger CM energies {y/s > 500 GeV) the fusion 
processes start to dominate over the ZH process. The ZZ-fusion process is 
suppressed by one order of magnitude compared to WW-fusion due to the 
ratio between neural and charged currents [6]. 

Both the signal and background events are produced at the centre of 
mass energy ^/s = 250 GeV, total integrated luminosity of 250 fb~^ and the 
Higgs boson is assumed to have mass 120 GeV. The choice of energy in this 
analysis maximizes the cross-section value for Higgs-strahlung. Standard 
Model events (mainly WW, ZZ and qq pairs) and Higgs decays to other 
fermions other than charm quarks are considered as backgrounds. All 0, 
2 and 4 fermion final states were generated using the Whizard Monte Carlo 
Event Generator [7]. PYTHIA |8] was used for the final state QED and QCD 
parton showering, fragmentation and decay to provide final-state observable 
particles. Photons from beamstrahlung and initial state radiation were also 
included in the simulations. 

For this study, event samples are created conforming to the expected 
ILC baseline parameters of ±80% electron and =f30% positron polarization. 
About 7M events were processed through the full detector simulation, with 
individual events weighted to reflect the statistical sampling. 

The detector response to generated events is simulated using the Geant4 
toolkit [HIITO], which provides the necessary classes to describe the geometry 
of the detector, the transport and interaction of particles with materials 
and fields. A thin layer of Linear Collider specific code, SLIC [llj, provides 
access to the Monte Carlo events, the detector geometry and the output of 
the detector hits. The detector parameters can be varied without having 
to rebuild the simulation executable binaries since the geometries are fully 
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Figure 2: Production cross-section for ZH process (solid lines) and all Hff 
processes (dotted lines) for a 120-200 GeV SM Higgs as a function of the 
centre-of-mass energy 

described at runtime. The output is in the standard LCIO format [12J so 
that detectors using other simulation packages can be analyzed and data 
generated using this system can be analyzed in other analysis frameworks. 

This analysis uses the Silicon Detector (sid02) model developed by the 
Silicon Detector (SiD) Design Study [13]. The tracking, electromagnetic 
calorimetry, highly segmented hadronic calorimetry using PEA (Particle Flow 
Algorithm) and vertex detector are all silicon based. The SiD also incorpo- 
rates a high field solenoid, iron flux return, and a muon identification system. 

The vertex detector (VD) is especially important for this analysis since 
it relies on flavour tagging based on the lifetime information of decaying 
particles. The VD comprises of a central barrel system with 5 silicon pixel 
layers and forward and backward disk regions each with 4 pixel disks. It 
is enclosed within a double walled carbon fibre support. The pixels are 20 
X 20/im^ and the sensor thickness is taken to be 75/im. Sensor cut width is 
8.68 mm in the innermost layer and 12.58 mm in all other layers. The sensor 
cut length is 125 mm for all layers. Carbon fibre cylinders support the sensors 
and the vertex detector is attached to the support tube with 25% material 
coverage support disks. For SiD, the charged track momentum resolution is 
expected to be better than a{l/pT) = 5 x 10~^{GeV/c)~^ for high momentum 
tracks and the impact parameter resolution will be better than ar<i) = (Jrz = 
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5 © 10/(psin^'^ 6*) [/iin]. The overall material budget for the VD per layer 
is 0.15% of the radiation length (Xq). Detailed geometry description of the 
detector is available at http://lcsim.org/detectors/7^sid02. 



3 Analysis Tools 

Identification of jets is an importand part of this analysis. The fragmen- 
tation products of hadronic systems are forced in N jets. The DURHAM 
algorithm [T3] is used in the analysis and provides y-cut values that deter- 
mine the separation between 1 and 2 jet cases and generally between N-1 
and N jet cases. 

To provide the most probable kinematic configuration of the event topol- 
ogy a kinematic fitter (Marlin Kinfit [15]) with four-momenta and mass con- 
straints is used. The fitter uses the method of Lagrange multipliers to deter- 
mine the most probable value for the jet four- momentum. 

In order to identify primary, secondary and tertiary vertices the topologi- 
cal vertex finder ZVTOP is used. The algorithm is part of a vertexing package 
developed by the LCFI collaboration [16]. It classifies events on the number 
of found vertices and combines eight optimized variables for each type of 
event in a neural network, which is then separately trained on samples of 
b-, c- and light quarks. The best discriminating variables are the corrected 
vertex mass, secondary vertex probability, impact parameter significance of 
the most significant track and the number of vertices in the event. There are 
nine networks used with eight inputs, one hidden layer with 14 neurons and 
one output neuron. 

Figure |3] shows the performance of the LCFI package optimized for the 
SiD detector. The plot shows the dependence of the purity on efficiency for 
the di-jet sample at y/s = 500 GeV for b-tagging, c-tagging and c-tagging 
with b-only background. The samples used for training of neural networks 
pass through the full SiD simulation and reconstruction. 

The open source Fast Artificial Neural Network (FANN) package p!7] is 
used to provide neural networks for the final event selection. The package is 
written in C programming language called from within a C-I--I- wrapper. 

4 Event Selection 

The analysis signature is dependent on the Z boson decay products (charged 
leptons, hadrons or neutrinos). The channels studied in this analysis are the 
neutrino mode (Z decaying to neutrinos) and the hadronic mode (Z decaying 
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Figure 3: Purity versus Efficiency for a 500 GeV di-jet sample for b-tagging 
c-tagging and c-tagging with b-only background. 

to hadrons). 

The selection of signal events is performed in three stages. The ffist step 
involves the classification of events into two channels using the number of 
leptons and visible energy in the event. Visible energy is defined as the sum 
of energies of all reconstructed particles in the event. Leptons are defined 
as reconstructed electrons or muons with momentum larger than 15 GeV. 
Figure H] shows the distributions of the visible energy and the number of 
leptons for the signal (ffiled histogram), Higgs background (dashed line) and 
SM background (solid line) before any selections. The neutrino mode is 
selected as events with no leptons and with visible energy in the 90 to 160 
GeV interval. The hadronic mode is selected as events with no leptons and 
the visible energy above 170 GeV. 

The second step is a cut based selection which reduces the backgrounds 
in the selected channel followed by a final stage using two neural network 
(NN) discriminants trained on signal and two types of background samples, 
SM and Higgs. The remaining events after the NN selection are used for the 
calculation of the cross sections and branching ratios. 

4.1 The neutrino channel 

In this channel all reconstructed particles are clustered into two jets which, 
for the signal, are assumed to come from the Higgs boson recoiling against 
two neutrinos from the Z boson decay. The discrimination between the signal 
and background with a different number of jets is achieved by selection of the 
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Figure 4: Visible energy (a) and the number of leptons per event (b) used for 
channel classification. Solid curves are SM background, dashed curves are 
inclusive Higgs sample and filled histograms are the signals. 

Vmin parameter. It corresponds to the minimum y-parameter for the two jet 
hypothesis in the Durham algorithm. The y-parameter is defined as follows: 

1. Given all pairs of particles i,j calculate 

2mm(E2^E|)(l-cos%) 

y^J = ^2 ' (1) 

where E^g is the sum of the visible energies of all particles before any 
recombination. 

2. If all yij > y-parameter stop calculation. The number of jets in the 
event is then defined as equal to the number of particles left. 

3. Otherwise recombine the pair with the smallest yij into a single particle. 

4. Finally, go back to step 1. 

Therefore, in this channel the ymin parameter is defined as the minimum 
value of the y-parameter that leads to a 2-jet event. Another powerful dis- 
criminant, the reconstructed invariant mass of two hadronic jets, is expected 
to be consistent with the Higgs mass. 

Figure [5] shows distributions of various variables used for pre-selection for 
the signal and backgrounds after classification. The kinematic variables used 
are: 

• Transverse momentum of jet, Pt'- The majority of SM background 
events are soft compared to signal events. The signal peak is expected 
to be at mi^/2. 
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• ^tracks PGf jet: Require more than 4 charged tracks to reject purely 
leptonic events 

• — log(ymm): This variable is used to exclude fully hadronic WW (and 
ZZ) events which are 4-jet events mis-identified as 2-jet events. The 
— log(?/mm) value is larger for the backgrounds compared to the signal. 

• Thrust: This is the maximum directed momentum given by 

T = max (2) 

" 2^\Pi\ 

i 

Signal events are less boosted and are more spherical than the back- 
ground events. 

• \cos{OthTust)\'- Signal events are produced more centrally in the detector 
while majority of the background processes have a strongly forward 
peaking angular distribution. 

• Angle between jets: A large fraction of background events have back- 
to-back jets but there also backgrounds with low angular separation 
due to low boost. Signal events are confined are confined to a range of 
angular values due to constraints coming from Z and H production. 

• Di-jet invariant mass: The signal di-jet mass is required to be between 
100 and 140 GeV to reject low mass hadronic systems from WW, ZZ 
and two-photon events 

• Highest reconstructed photon energy: The signal events rarely have 
hard photons and this helps reject both highly energetic initial state 
radiation (ISR) photons and hard photons from di-jet events which 
occur for the background. No photon isolation is required. 

The full list of selections is given in Table [TJ 



The background includes all SM processes. The most important are 2- 
fermion events, ZZ pairs decaying to neutrinos and hadrons and WW pairs 
where one W decays hadronically and the other W decays into a neutrino 
and a lepton which escapes undetected. Table [2] shows the number of events 
before and after pre-selection cuts. 

The remaining events are categorized using the neural networks imple- 
mented in FANN [T7] . The NN input variables include all the variables stated 
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Figure 5: Distribution of pre-selection kinematic variables in the neutrino 
channel. Solid curves are SM background, dashed curves are inclusive Higgs 
sample and filled histograms is the signals. 
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Cut 






selection 


value 


1. 


20 


< 


Pt of jet 


< 90 GeV 


2. 






number of charged tracks per jet 


> 4 


3. 






- log(2/mjn) 


< 0.8 


4. 






thrust 


< 0.95 


5. 






COs{6thrust) 


< 0.98 


6. 


100° 


< 


angle between jets 


< 170° 


7. 


100 GeV 


< 


di-jet invariant mass 


< 140 GeV 


8. 






Highest reconstructed photon energy 


< 10 GeV 



Table 1: Selections for the neutrino channel 



Cuts 


SM background 


Higgs background 


Signal 


(i) Before Classification 


9275594683 


17475 


652 


(0) After Classification 


45936973 


14294 


637 


(1) 


18374789 


13783 


619 


(2) 


17123140 


13729 


618 


(3) 


6849256 


13416 


609 


(4) 


685329 


12179 


560 


(5) 


627113 


11945 


550 


(6) 


576422 


10226 


525 


(7) 


203292 


10088 


519 


(8) 


109057 


9902 


506 



Table 2: Number of cc events before channel classification, after channel 
classification and after pre-selection cuts in the neutrino mode. 



above with addition of the LCFI flavour tag outputs for both jets. Figure 
shows the distribution of the three possible LCFI flavour tags, 'b-tag', 'c-tag' 
and 'c-tag with b background only' for the leading jet. The first NN is trained 
to distinguish the SM background from the inclusive Higgs sample and to 
produce the '^'^SM-Higgs output; and, the second NN to distinguish the sig- 
nal from the Higgs background sample and to produce the ^^mggs-signai 
output. The outputs of the trained NNs are shown in Figure [71 

Figure [8] shows scatter plots of ^^Higgs-signai versus NNsM-Higgs for the 
signal and backgrounds, and also shows the number of events in different 
regions of Higgs-signai and '^'^sM-Higgs parameter space. The final event 
sample is required to have Higgs-signai > 0.2 and ^'^SM-Higgs > 0.3. This 
cut is optimised by using different values of ^^mggs-signai and ^'^SM-Higgs 
and choosing a point at which the signal cross-section error is minimised (or 
S/ \/B ratio is maximal). The optimization is shown Figure |9l After the final 
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Figure 6: Neutrino channel: (a) b-tag; (b) c-tag; (c) c-tag with b background 
only;. Solid curves are SM background, dashed curves are Higgs background 
sample and filled histograms are the signal. 

NN signal weighted | NN2 signal weighted 1 




(a) (b) 

Figure 7: Neutrino channel: (a) NNgM-Higgs output; (b) NN^iggs-signai out- 
put. Solid curves are SM background, dashed curves are Higgs sample and 
filled histograms are the signals. 
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selections the sample includes 178 signal events with SM background of 140 
events and Higgs background of 109 events. The signal efficiency is about 
28% and the purity is 42%. The efficiency is calculated as the ratio between 
the final number of signal events and the number of signal events after the 
classification but before the pre-selection cuts. 
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Figure 8: ^^Higgs-signai versus 1^1^ SM -Higgs foT (a) Signal, (b) Standard 
Model background, (c) Higgs background and (d) number of events for dif- 
ferent NN regions 



4.2 The hadronic channel 

In the hadronic channel, events are forced to have four reconstructed jets 
and as in the two-jet case the i/min variable is chosen to differentiate from 
backgrounds with different numbers of jets. In this channel the ymin param- 
eter is the minimum value of the y-parameter that leads to a 4-jet event 
configuration. For the signal, two of the jets should have di-jet invariant 
mass consistent with the Higgs boson and the other two having di-jet mass 
consistent with the Z boson. 
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(a) 



Figure 9: Neural Network Cut optimization for cc in the neutrino channel. 

Jet pairing is performed before kinematic fitting. For the 4-jet events we 
have 6 possible pairings of the jets and 3 possible associations of the 4 jets 
to the Z and H bosons. For the 6 possible pairings we calculate the invariant 
mass of each pair and compare with the masses of bosons. For each event we 
calculate 



The pairing that minimises d is chosen and we find on average that jets 1 
and 3 are associated to the Higgs and jets 2 and 4 are associated to the Z 
boson, where the ranking is performed according to jet energy (jet 1 has the 
highest energy). 

Kinematic fitting |T5] is performed to improve identification of two jets 
from Higgs boson and two jets from Z-boson in order to reduce the combina- 
torial background. The kinematic fit minimizes a global chi square which is 
obtained from the difference between fitted and measured parameters of the 
kinematic variables, which are weighted by their covariance matrices. The 
following constraints are used for the analysis: 




(3) 



3 



i=0 



• sum of E 



250 GeV 
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• invariant mass difference of two jet pairs is equal to difference between 
Higgs and Z boson masses 

The main backgrounds for this channel are WW and ZZ pairs where the 
all bosons decay to hadrons. Figure dO] shows the Higgs and Z invariant 
masses before and after kinematic fitting. 



H mass before and after fitting~l Z mass before and after fitting~| 




(a) (b) 



Figure 10: Higgs (a) and Z (b) invariant masses before (dashed line) and 
after (solid line) kinematic fitting. 

For the preselection, we use kinematic cuts similar to the cuts for the 
neutrino channel. Figure [11] shows distributions of variables used for pre- 
selection for the signal and backgrounds after classification. The variables 
are as defined in Section 4.1. 

Selections used for further reduction of the background are presented in 
Table O For the cuts, the jets were ordered in energy and the invariant mass 
and angle cuts were employed after the kinematic fitting. Table H] shows the 
number of events before and after pre-selection cuts. 

For a 250 GeV centre-of-mass scenario, the energy of the photon resulting 
from radiative returns to the Z-peak is expected to be at around 108.6 GeV. 
Figure fTU i) shows the distribution of the SM background reconstructed pho- 
ton energy after channel classification but before any kinematic selections. 
There is a peak seen at about 125 GeV which does not correspond to the 
radiative return to Z-peak. The main contributions to this peak are events 
where a colliding e"*" or e~ loses most of its energy by ISR. Figure [T^ a) 
shows the energy of all SM background photons at the generator level for 
the hadronic channel after the channel classification. This distribution has 
two peaks at about 108 GeV and 125 GeV, which confirms that these pho- 
tons originate from both Z radiative returns and direct ISR. The first peak 
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Figure 11: Distribution of pre-selection kinematic variables. Solid curves 
are SM background, dashed curves are Higgs background sample and filled 
histograms are the signal. See text for discription of variables. 
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corresponds to the radiative return to the Z-peak and the second peak cor- 
responds to the peak observed in Figure fTTT i). Figure [TWb) shows the SM 
reconstructed photon energy after channel classification but here all kine- 
matic selections (1-8) are applied. Here we notice that we now have a peak 
at 108 GeV and no peak at around 125 GeV. We conclude that the previous 
selections are successful in removing direct ISR photons in the other peak. 



MC Photon without parents energy after classification | Photon energy after all culs before Ephol cut I 




Figure 12: (a) Energy of MC photon for the Standard Model background 
after channel classification; (b) Energy of reconstructed photon for the Stan- 
dard Model background after channel classification and kinematic selection 
cuts. 



Cuts 






selection 


value 


1. 






number of charged tracks per jet 


> 4 


2. 






- log(l/mm) 


< 2.7 


3. 






thrust 


< 0.95 


4. 






COs{6thrust) 


< 0.96 


5. 


105° 


< 


angle between jet 1 and 3 


< 165° 


6. 


70° 


< 


angle between jet 2 and 4 


< 160° 


7. 


110 GeV 


< 


invariant mass of Higgs candidate after fit 


< 140 GeV 


8. 


80 GeV 


< 


invariant mass of Z candidate after fit 


< 110 GeV 


9. 






Highest reconstructed photon energy 


< 10 GeV 



Table 3: Selections for the four-jet analysis. 



Similarly to the neutrino channel the variables above and also the three 
fiavour tag outputs for all jets, are used in a neural network based selection 
employing the FANN package. Figure [13] shows the distribution of the three 
LCFI fiavour tags, 'b-tag', 'c-tag' and 'c-tag with b background only' for 
the first and second leading jets. The first NN is trained to distinguish the 
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Cuts 


SiVl background 


TT" "U 1 J 

Higgs background 


* 1 

Signal 


(i)Beiore Llassiiication 


9/75594o83 


51513 


1931 


(UjAiter Ulassincation 


39398366 


41016 


1837 


(1) 


1 opm Tlx o 

106U1753 


on 1 C 

30125 


1 1 /1 
1143 


(2) 


13921271 


29478 


1101 


(3) 


8737017 


27065 


1047 


(4) 


7943851 


26322 


1017 


(5) 


5871237 


26001 


979 


(6) 


4898312 


25687 


978 


(7) 


1917231 


22533 


966 


(8) 


1561432 


21877 


963 


(9) 


967312 


15687 


947 



Table 4: Number of cc events before channel classification, after channel 
classification and after pre-selection cuts in the hadronic mode. 



SM background from the inclusive Higgs sample producing the '^^sM-mggs 
output, and, the second NN is trained to distinguish the signal from the 
Higgs background sample producing the "^"^uiggs-signai output. The outputs 
of the trained NNs are shown in Figure [TH 

Figure [I5] shows scatter plots of ^^mggs-signai versus ^^sM-mggs for 
the signal and backgrounds, and also shows the number of events in different 
regions of ^^Higgs-signai and '^^^'^sM -Higgs- The final signal sample is required 
to have "^^uiggs-signai > 0.2 and '^'^sM-Higgs > 0.3. Optimisation of the NN 
cuts used the same method as in the neutrino channel. The optimization 
is summarized in Figure [161 After the final selections the total number of 
events sample includes 407 signal events with SM background of 673 events 
and Higgs background of 213 events. The signal efficiency (calculated as in 
the neutrino channel) is about 22% and the purity is 32%. 

5 Branching ratio calculation and results 

The branching ratio of the Higgs boson decay to charm quarks was calculated 
using events that passed the final neural network selection. The calculation 
was done by normalising the signal cross section to the inclusive Higgs cross 
section, azH = 209±9.8 fb, as determined in an independent recoil mass 
analysis performed for the SiD Letter of Intent [13] . The branching ratio is 
then given by 

BRiH ^ //) = ^ (4) 
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I total tagged bs jet 1 after cuts | 



I total tagged cs jet 1 after cuts 




Figure 13: Hadronic channel: (a) first jet b-tag; (b) first jet c-tag; (c) first jet 
c-tag with b background only (d) second jet b-tag; (e) second jet c-tag; (f) 
second jet c-tag with b background only. Solid curves are SM background, 
dashed curves are Higgs background sample and filled histograms are the 
signal. 
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Figure 14: Hadronic channel: (a) l^'NsM-Higgs output; (b) N^mggs-signai 
output. Solid curves are SM background, dashed curves are inclusive Higgs 
sample and filled histograms are the signals. 
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Figure 15: ^'^Higgs-signai versus NNsM-Higgs for (a) Signal, (b) Standard 
Model background, (c) Higgs background and (d) number of events for dif- 
ferent NN regions. 
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(a) 



Figure 16: Neural Network Cut optimization for cc in the hadronic channel. 



where / represents the daughter decay products from the Higgs. The relative 
accuracy of the the branching ratio takes into account both the relative signal 
cross section uncertainty and the relative Higgs-strahlung uncertainty given 
as 

ABR U^V(Aa^' 



BR 

with the relative signal cross section uncertainty calculated by 



Aanff V signal + background 
(Jh f f signal 

and the cross-section is calculated as follows 



(6) 



^Hff = f (7) 

where N is the number of signal events after all selections, e is the efficiency 
of signal selection and L is the total integrated luminosity. 

The uncertainty in the efficiency is considered negligible, relying on sim- 
ulations to determine it with sufficient precision. The systematic effects or 
contributions of the luminosity uncertainty were not considered in this anal- 
ysis. 
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The weighted average of the signal cross section and its uncertainty are 
calculated using cross section and relative uncertainty values obtained from 
the neutrino and hadronic channels. The weighted average cross section is 
given by 

_ x{6yY + y{6x)^ 

^average- (5^)2 + (5^)2 ' 

where x and y are the cross sections in the neutrino and hadronics channels 
respectively, and 5x and 5y are the cross section uncertainties in the neutrino 
and hadronic channels respectively assuming that the two channels are sta- 
tistically independent. The uncertainty of the average cross section is then 
calculated as 

5x * 5y 



^{5xY + {5yf 



(9) 



where 5z is the uncertainty of the weighted average cross section. 

A summary of the results obtained for both the neutrino and hadronic 
channels is given in Tables [5] The relative accuracy of the branching ratio in 
both channels is dominated by the precision of determination of the signal 
cross section since the uncertainty of the total Higgs-strahlung cross section is 
small (4.7%). The larger value of the BR uncertainty in the neutrino channel 
is largely due to the lower signal statistic since BR(Z — )■ up) is 20.0% ,and 
BR(Z ^ qg) is 69.91%. 





Neutrino 


Hadronic 


Combined 


Signal events 


178 


407 




SM background events 


140 


673 




Higgs background events 


109 


213 




Signal efficiency % 


27.943±0.018 


22.156±0.010 




Signal aHcc 


6.8±0.8 fb 


6.9±0.6 fb 


6.86±0.48 fb 


Relative uncertainty on aucc 


11.6% 


8.8% 


7.0% 



Table 5: Measurement results of H— )■ cc branching ratio. 



Using Equation 5 the uncertainty on the branching ratio is found to be 
.4%. 



6 Conclusion 

The uncertainty of the measurement of the decay branching ratio to charm 
quarks for a neutral SM Higgs boson of mass 120 GeV has been studied at a 
centre-of-mass of energy of ^/i = 250 GeV and a total integrated luminosity 
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of 250 y fb ^. The analysis is based on full detector simulation and realistic 

event reconstruction. The relative uncertainties obtained are comparable to 
values obtained in some of the previous studies. A good performance of 
flavour tagging and the use of neural networks in event selection are critical 
in obtaining these results. In this study we find the uncertainty on the cross 
section to be 11.6% in the neutrino channel and 8.8% in the hadronic channel 
giving a combined uncertainty on the branching ratio of 8.4%. 
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